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Abstract. We present the analysis of the source counts in the XMM-COSMOS 
survey using data of the first year of XMM- Newton observations. The survey 
covers ~2 deg 2 within the region of sky bounded by 9 h 57.5 m < R.A. < 10' l 03.5 m ; 
l d 27.5 m < DEC < 2 d 57.5 m with a total net integration time of 504 ks. Using a 
maximum likelihood algorithm we detected a total of 1390 sources at least in one 
band. Using Monte Carlo simulations to estimate the sky coverage we produced 
the logN-logS relations. These relations have been then derived in the 0.5-2 keV, 
2-10 keV and 5-10 keV energy bands, down to flux limits of 7.2 x 10 -16 erg cm~ 2 
s _1 , 4.0xl0 -15 erg cm -2 s _1 and 9.7xl0~ 15 erg cm" 2 s _1 , respec tively. These 
relations have been compared to previous X-ray survey and to the most recent 
X-ray background model finding an excellent agreement. The slightly different 
normalizations observed in t he source counts of COSMOS and previous surveys 
can be largely explained as a co mbination of low counting statistics and cosmic 
variance introduced by the large scale structure. 



1. Introduction 



A solid knowledge of the X-ray source counts is fundamental to fully under 
stand the AGNs population of the X-ray background (XRB). At present the 
two d eepest X-ray surveys, the Chandra Deep F ield North (CDFN; iBauer et al 



20041 ) and Chandra Deep Field South fCDFS: iGiacconi et al.ll200ll h have ex- 



tended the sensitivity by about two orders of magnitude in all bands with re- 
spect to previous surveys ( Hasinger et al.|[l993l ; lUeda et al.|[l999l ; iGiommi et ail 
2000), detecting a large number of faint X-ray sources. However, deep pencil 



beam surveys are limited by the area which can be covered to very faint fluxes 
(typically of the order of 0.1 deg 2 ) and suffer from significant field to field vari- 
To ovecome this problem, shallower surveys over larger areas have been 



ance. 



underta ken in the last few years with both Chandra (e.g. the 9 deg^ Bootes 
surve y ( Murray et al.l [20051 ) . the Extended Groth strip EGS (Nandra et al 



200511. the Extend ed Chandra Deep F ield South E-CDFS. dLehmer et a"l l2005l: 

4 



Virani et all 120061 ) and th e Champ ([Green et al.1 120041: [Kim et alJ l2004n1 and 
XMM-Newton (e.g. the HELLAS2XMM survey ([Fiore et alJ l2003l). the XMM- 
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Newt on BSS ([Delia Ceca et al.ll2004h and the ELAIS SI survey (|Puccetti et al 
2006) ). I n this context the X MM-Newton wide field surve y in the COS- 



MOS field (jScoville et alJl200d ). hereinafter XMM-COSMOS (jHasinger et al 
2006), has been conceived and designed to maximize the sensitivity and sur- 



vey area product, and is expected to provide a major step forward toward a 
complete characterization of the physical properties of X-ray emitting Super 
Massive Black Holes (SMBHs). In this work we concentr ate on the first year o f 
XMM-Newton observations (A03) of the COSMOS field (jHasinger et al.ll2006l ). 
For more details see Cappelluti et al. (2007). 



2. The X-ray logN-logS 

The source detection was performed using a maximum likelihood algorithm pro- 
vided with the XMM-Newton Standard Analysis System (SAS) in the 0.5-2 
keV, 2-10 keV and 5-10 keV bands. A total of A total of 1281, 724 and 186 
point-like sources were detected in the three bands down to limiting fluxes of 
7.2xl0~ 16 erg cm~ 2 s -1 , 4.7xl0~ 15 erg cm -2 s _1 and 9.7xl0~ 15 erg cm -2 
s _1 , respectively. Using Monte Carlo simulation to determine the sky cover- 
age, source number counts can be easily computed using the following equation: 
N(> S) = X)j=i T[~deg~ 2 , where N$ is the total number of detected sources 
in the field with fluxes greater than S and f2j is the sky coverage associated 
with the flux of the \ th source. The cumulative number counts, normalized to 
the Euclidean slope (multiplied by S ), are shown in Figure 1, for the 0.5-2 
keV, 2-10 keV and 5-10 keV energy ranges, respectively. We performed a bro- 
ken power-law maximum likelihood fit to the unbinned differential counts in the 
0.5-2 keV and 2-10 keV energy bands. In the 0.5-2 keV energy band the best 
fit parameters are ai=2.60lg;lf| ) a 2 =l-65±0.05, S 6 =1.55±g;i xl0~ 14 erg cm" 2 
s _1 and ^4=123. Translating this value of the normalization to that for the 
cumulative distribution at 2xl0~ 15 erg cm~ 2 s _1 , which is usually used in the 
literature for Chandra surveys, we obtain A15 ~450 which is fully consistent with 
most of previous works where a fit result is presented. In the 2-10 keV band 
the best fit parameters are «i=2.43±0.10, a 2 =l-59±0.33, S b =1.02+^ xl0~ 14 
erg cm -2 s _1 and ^4=266. The latest value translates into ^15 ~1250. Also in 
this band, our results are in agreement with previous surveys within la. In the 
5-10 keV energy bands, where the differential counts do not show any evidence 
for a break in the sampled flux range, we assumed a single power-law model 
for which the best fit parameters are found to be ^4=102 and ai=2.36±0.1. 
In the 0.5-2 keV band we measure a contribution of the sources to the XRB 
which corresponds to a normalization at 1 keV of 7.2±1.1 keV cm -2 s" 1 keV -1 . 
The corresponding values in the 2-10 keV and 5-10 keV bands are 4.7±1.1 keV 
cm" 2 s- 1 keV- 1 and 2.6 ±1.4 keV cm" 2 s" 1 keV" 1 . Therefore XMM-COSMOS 
resolves by itself ~65%, ~40% and ~22% of the XRB into discrete sources in 
the 0.5-2 keV, 2-10 keV and 5-10 keV energy bands, respectively. In Figure 
1 we compared our logN-logS to t hose predicted by the recent XRB model of 
Gilli. Comastri fc Hasinger I (|2006h . 



The amplitude of source counts distributions varies significa ntly among dif- 
ferent surveys (see e.g. lYang et "al1l2003l : ICappelluti et al. 1120051 . and references 



therein). This "sample variance", can be explained and predicted as a com- 
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Figure 1. The X-ray logN-logS compared to the survey listed in the 
labels in the 0.5-2 keV, 2-10 keV and 5-10 keV energy bands, respec- 
tively Top Left, Top Right and Bo ttom Left Panels. The b lack — 
dashed line represents the prediction of lGilli. Comastri &: Hasinger I ((2006) . 
Bottom Right Panel : The counts in cell fluctuations within the XMM- 
COSMOS field on different angular scales. The dashed lines represent the 
la expected fluctuation. 



bination of Poissonian variation s and effects due to the clustering of sources 
( Peebles 19801 ; Yang et al. 20031 ). In order to determine whether the differences 
observed in the source counts of different surveys could arise from the clustering 
of X-ray sources, we estimated the amplitude of the fluctuations from our data, 
by producing subsamples of our survey with areas comparable to those of. e.g., 
Chandra surveys. 

The XMM-COSMOS field and the Monte Carlo sample fields of Section 4 were 
divided in 4,9,16 and 25 square boxes. Making use of the 0.5-2 keV energy band 
data, we computed for each subfield, the ratio of the number of real sources to 
the number of random source. Both the random and the real sample were cut 
to a flux limit of 5xl0 -15 erg cm -2 s _1 . In the lower right panel of Figure 1 
we plot the ratio of the data to the random sample as a function of the size of 
the cells under investigation. The predicted fractional standard deviations are 
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0.13, 0.19, 0.23 and 0.28 on scales of 0.44 deg 2 , 0.19 deg 2 , 0.11 deg 2 and 0.07 
deg 2 , respectively. The measure fractional standard deviations are 0.09, 0.20, 
0.21 and 0.24 on the same scales, respectively. The ratio of the dusting to the 

Poissonian variance is expected to scale as o c ija v oc A/" ' 5 ^ 7 1 ^ //2 a^ 3 ~ 7 ^ 2 . We 
therefore conclude that at this scales and at fluxes sampled here the major con- 
tribution to source counts fluctuations is due to Poissonian noise. This analysis 
is at least qualitatively consistent with Figure 2, which shows a significantly 
larger dispersion in the data from different surveys in the hard band than in 
the soft band. Moreover, the results here discussed are also consistent with the 
observed fluctuations in the deep Chandra fields (see, for example, Bauer et al. 
2004). Large area, moderately deep surveys like XMM-COSMOS are needed 
to overcome the problem of low counting statistics, typical of deep pencil beam 
surveys, and, at the same time, to provide a robust estimate of the effect of large 
scale structure on observed source counts. 
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